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Premicellar aggregation of amphiphilic molecules: 
Aggregate lifetime and polydispersity 
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A recently introduced thermodynamic model of amphiphilic molecules in solution has yielded, 
under certain realistic conditions, a significant presence of metastable aggregates well below the 
critical micelle concentration — a phenomenon that has been reported also experimentally. The 
theory is extended in two directions pertaining to the experimental and technological relevance of 
ON \ such premicellar aggregates, (a) Combining the thermodynamic model with reaction rate theory, 

we calculate the lifetime of the metastable aggregates, (b) Aggregation number fluctuations are 
examined. We demonstrate that, over most of the metastable concentration range, the premicellar 
aggregates should have macroscopic lifetimes and small polydispersity. 
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I. INTRODUCTION 

The natural and technological applications of self-assembled amphiphilic structures (micelles) in aqueous solution 
are vast [H, 0]. According to the common view, supported by numerous macroscopic experiments (e.g., conductivity 
and surface-tension measurements) and widely accepted theories 0, Ej| , amphiphilic molecules form aggregates above 
a well defined critical micelle concentration (cmc). During the years, however, there have been several experimental 
indications [1, H, S S @] ; as well as theoretical ones @ , for the appearance of aggregates at concentrations well below 
the cmc — aphenomenon referred to as premicellar aggregation. In particular, a fluorescence correlation spectroscopy 
experiment [6| seems to have provided direct observation of premicellar aggregates at concentrations four times lower 
O \ than the macroscopically determined cmc. 

O . Recently we have presented a two-state (monomer-aggregate) thermodynamic model for amphiphilic aggregation 
which, alongside its simplicity, allows the study of metastable micelles of variable size The analysis has yielded 
£NJ ■ a sequence of three well separated concentrations: ci, where a metastable aggregated state appears but is not signifi- 
' cantly occupied; c%, above which an appreciable amount of metastable aggregates forms; and C3, where the aggregated 
state becomes stable. The cmc as commonly measured in macroscopic experiments has been shown to correspond 
to C3. Thus, appreciable premicellar aggregation may occur in the concentration range between Ci and C3. We have 
shown that, so long as the micelles are not too large, the extent of premicellar aggregation is much larger than what 
would be expected from mere finite-size effects. This somewhat surprising effect stems from the variability of the "ex- 
cited state", i.e., the freedom of the micelles to select their sizes, and the small free-energy difference between the pure 
monomeric state and the metastable one, which contains mostly monomers and a low concentration of aggregates. In 
addition, we have found that the premicellar regime is characterized by a weak concentration dependence of micelle 
size. Thus, the premicellar aggregates have roughly the same size as those observed above the cmc, in agreement with 
the experiment of Ref. |g|. 

The analysis in Ref. [lOj, which is purely thermodynamic, does not fully account for possible kinetic barriers for 
premicellar aggregation. It has been assumed, and will be assumed below, that the solution has fully equilibrated. 
(In cases where high nucleation barriers exist, they might be overcome in practice, e.g., through heterogeneous 
nuclcation.) Two key issues remain open, however. First, while the metastable premicellar state may be appreciably 
occupied at equilibrium, the aggregates might be short-lived. Second, although the mean size of the premicellar 
aggregates is similar to that of the micelles above the cmc, the size distribution in the former case might be much 
broader. Evidently, these issues of lifetime and polydispersity could jeopardize the experimental and technological 
relevance of premicellar aggregation. 

Micellization dynamics (above the cmc) were thoroughly studied in previous works, both experimentally (see Ref. 
[TH and references therein) and theoretically (see, e.g., Refs. [H, EH, 0, EE EE E3) EH EE H3] ) • Two disparate time 
scales are involved in the dynamics, corresponding to the exchange of individual monomers between the micelle and 
the solution and the much slower process of micelle formation and breakup. Being interested in aggregate stability, we 
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focus here on the latter. We use the free energy landscape, as obtained from the thermodynamic model [lOj, within 
Kramers' rate theory [2l], to study the lifetime of premicellar aggregates. The second extension of the theory is 
an examination of aggregate size fluctuations in the premicellar regime. 

In Sec. [II] we briefly review the thermodynamic model of Ref. [101 ] and then extend it to study lifetime and poly- 
dispersity. Representative numerical results are presented in Sec. IIIII In Sec. [IV] we discuss the results and their 
implications for the observation of premicellar aggregates. 



II. MODEL 



A. Free energy 

Our starting point is the two-state thermodynamic model of micellization presented in detail in Ref. [nj. The 
solution, containing a volume fraction <fi of amphiphiles, is assumed to consist of two species: monomers with volume 
fraction </>i, and aggregates of m molecules with volume fraction (</3 — <f>i). The volume fraction of water is (1 — 0). 
Both 4>i and m are treated as degrees of freedom, i.e., the system can select the number of aggregates as well as their 
size, while the total volume fraction cj> is the control parameter. The two-state approximation restricts the validity of 
the entire approach to compact (spherical) micelles, whose size distribution is relatively narrow [2J. The free energy 
of the solution contains a mixing-entropy contribution and an interaction term. The former is calculated using a 
coarse-grained (Flory-Huggins) lattice scheme, where a water molecule occupies a single lattice site (of volume a 3 ), 
and each amphiphile occupies n sites. The latter term, containing all other contributions to the free energy of transfer 
of a monomer from the solution to an aggregate of size m, is represented by a single phenomenological function, u(m). 
The resulting free-energy density (per lattice site) is 

F B (0i,m, <t>) = — In 0i + — — pn(0 - <pi) ~ mu(m)} + (1 - 4>) ln(l - 0). (1) 
n nm 

(All energies in this paper are expressed in units of the thermal energy k^T.) The specific form of u{m) is not crucial 
for the analysis; it should merely have a maximum at a finite value of m to ensure the formation of finite aggregates 
(rather than a macroscopic phase) upon increasing cj>. For the sake of numerical examples we shall use the following 
function 0, El: 

u(m) =Uq — cm -1 / 3 - kto 2 / 3 . (2) 

(The physical origins of the terms appearing in Eq. (|2j) , as well as the rather limited range of relevant values for the 
parameters uq, ct, and k, are discussed in Ref. [To|.) 

Equation ([TJ) defines a free-energy landscape over a two-dimensional space of macrostates (<pi , m) . Along the <\>\ 
axis F s is always convex, i.e., it has a single minimum at 0™ m (m, 4>) for all values of m and for any <f> fl(]j |. Along the 
m axis the free energy becomes nonconvex above a certain volume fraction, <f> > c%, with two minima at m = 1 and 
m mm , and a maximum in-between, at m = m max . In the premicellar regime of interest, c\ < C2 < <f> < c 3j the free 
energy has a global minimum still at the pure monomeric state, [0™ m (m = l),m = 1], as well as a local minimum 
at the metastable aggregated state, [</f™ m (m mm ), 77i mm ] (containing mostly monomers and a a low concentration of 
aggregates). The two minima are separated by the saddle point [0™ m (m max ), m max ], which poses a kinetic barrier for 
the disintegration of the metastable aggregated state into the stable monomeric one. 

The following analysis relies on two basic assumptions. First, we assume that overcoming the barrier at the saddle 
point [</>f lm (m max ), m max ] is the rate-limiting process in aggregate dissociation, whereas diffusion is much faster. 
Hence, the dynamics depend on m alone, advancing at all times t along the path [^j" un (m(t)),m(t)]. The second 
assumption arises from the necessity to relate our coarse-grained model with single- aggregate properties. Since the 
model [e.g., Eq. |T])] does not explicitly consider single aggregates but rather macrostates containing both monomers 
and aggregates, we shall consider, instead, a fictitious subsystem, of volume V\, which on average contains a single 
aggregate of size m mm . The volume of the aggregate itself is na 3 m mm , and the volume fraction of aggregates is 4> — 4>i. 
Hence, the subsystem volume is 

na 3 m mi "(0) 

1W 0-0f n (m min (0))' 1 ' 

Since <f> — <pi is very small, V\ is far from being microscopic, and we may apply our coarse-grained description to the 
subsystem, writing its free energy as 

F(tf, 1 ,m,4>) = ^-F s , (4) 
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where F s is given by Eq. fl}. Thus, the dissociation of a single premicellar aggregate is treated as the transition of a 
mesoscopic subsystem from a metastable state, containing monomers and (on average) one aggregate, to the stable, 
purely monomeric state. For brevity the free energy of the subsystem along the dissociation path [4>™ in (m(t)) , m(t)} 
is hereafter referred to as F(m). 



B. Aggregate lifetime 

We follow the lines of Kramers' theory [2l|, while adapting it to the case of premicellar aggregates. The main 
assumptions of this approach are as follows, (i) The energy barrier between the two states is sufficiently high, leading 
to separation of time scales between the fast monomer exchange process and the much slower aggregate dissociation, 
(ii) The free energy of the final (monomeric) state is much lower than that of the initial (aggregated) one, ensuring 
a practically unidirectional probability current from the aggregated to the monomeric state. The first assumption 
breaks down when <fi is too small, i.e., as it gets too close to c\\ in the examples of Sec. IIIII it becomes invalid already 
for (f> ~ C2- The second assumption fails when cj) gets close to C3. Thus, the following calculation is strictly valid only 
for C2 C < C3. (The behavior outside this domain of validity will be commented on separately in Sec. II VI ) In 
addition, we assume that the aggregation number is large, m > 1, so that the discrete changes in m can be replaced 
to a good approximation by continuous, infinitesimal ones. 

We begin with the master equation for the probability density function, /(m, i), of finding the subsystem around 
the state [4>™ in (m), m] at time t, 



df(m,t) 
dt 



dkW{m - k, k)f(m -k,t)- / dkW(m, k)f{m, t)dk, 



(5) 



where W(m, k) is the transition probability per unit time for the aggregation number to change from mtoml k. 
Assuming that large jumps in aggregation number are improbable, we expand the first integral in Eq. ([5]) to second 
order in small k and get the Fokker-Planck equation, 



^ = - J^, j(m,t) = A(m)f(m,t) - ^-[D(m)f(m,t)}. 



(6) 



The first term in the probability current density j describes a drift along the aggregation-number axis, with velocity 
A(m) = J dkkWfa, k). The second term represents diffusion along that axis, with a diffusion coefficient 



D(m) = ~J dkk 2 W{m,k). 



Demanding that / reduce at equilibrium (i.e., when j = 0) to the Boltzmann distribution, f oq (m) 
gets from Eq. ^ a generalized Einstein relation between A and D, 

A(m) = -D(m)F'(m) + D'(m), 

where a prime denotes a derivative with respect to m. Substituting this relation back in Eq. 
probability current density as 



-F(r. 



(7) 



one 



j = -D{m)e- F ^ m \f{m,t)e F{m) ]'. 



(8) 

we rewrite the 
(9) 



Thanks to the assumed high free-energy barrier, and the resulting separation of time scales, steady state can be 
assumed practically throughout the entire dissociation process. Thus, df /dt — dj/dm — 0, i.e., j — j ss independent 
of m. Equation (J5|) can then be integrated over m, 



dm- 



D 



-fe 1 



(10) 



The second assumption, of a large free-energy difference between the two states, implies that the right-hand side 
(rhs) of Eq. (fT0|) is dominated by its value at to" 1111 . In addition, we assume that the subsystem is still mostly in the 
aggregated state near m min , at quasi-equilibrium, and, hence, f(m,t) ~ e~ F ^ m \ Expanding about m mm we obtain 
for the normalized probability density, 



f(m) ^ [F"(m min )/(27r)] 



1/2 p -5i r "(m min )(m-m mta ) 2 



(11) 
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The rhs of Eq. (TQ) is given, therefore, by -[F"(TO min )/(27r)] 1 / 2 e F ( mm " 1 ). 

Treating the left-hand side (lhs) of Eq. (fTO]) requires an estimate for the aggregation- number diffusion coefficient, 
D(m). We use the definition of this coefficient, Eq. (0, together with Langer's formula for the transition probability 

W(m,k) ~ T -i e -fe 2 /(2A) e -i[F( m +fc)-f( m )] ; (12) 

where To is a molecular time scale, and A is used to suppress large jumps in the aggregation number. Assuming that 
jumps much larger than unity are improbable, we set A = f . We then expand F(m + k) — F(m) in Eq. (fT2]) to second 
order in k, normalize the transition probability, and substitute it in Eq. ([7]) to obtain 



1 4 + F' 2 + 2F" 
2^ (2 + F") 2 



D{m) = — To , • (13) 



Analysis of Eqs. ([TJ) — (I3J) and (TT5)) shows that for realistic aggregation numbers, ra>l, one has \F\ ^> | ln£>|. Hence, 
the integral on the lhs of Eq. (jTUJ) is dominated by a small region around the maximum of F . We expand F(m) about 
m max , integrate, and get for the lhs of Eq. iJTOJ), j^[2n/\F"(m m ™)\] 1 / 2 e F ( mm ^ /D(m max ). 
Substituting all these results in Eq. (fT0|) . we finally obtain for the micelle lifetime (25l |. 



r D = bssT 1 = ^ ^e Fb , (14) 

(1 + F"(m max )/2) \F "(m min )F"(m max )| 1 

where Ft, = F (m max ) — F (m min ) is the height of the free-energy barrier between the aggregated and monomeric states. 
Equation (fl4|) . combined with Eqs. dU) — (SJ) , yields the aggregate lifetime in the metastable, premicellar regime. 



C. Polydispersity 

For a given amphiphilc volume fraction in the premicellar regime, C2 < 4> < C3 , the aggregation number of the 
metastable aggregates, m min (0), is given by the local minimum of F s of Eq. (fTJ) 10]. To examine the polydispersity 
of the aggregates we should calculate the fluctuations of m around m mm for a single aggregate. As explained in Sec. 
Ill A[ within our coarse-grained framework we calculate, instead, the fluctuations of m in a subsystem of volume V\. 
The distribution of m in that subsystem is given, for m close to m mm , by Eq. (fTTj) . Thus, we readily get for the 
mean-square size fluctuation, 

(5m 2 ) = 1/F"(m nlin ). (15) 

The relative width of the size distribution, 

(5m 2 ) 1 / 2 1 
w = = (16) 

(m) m min [F"(m min )] 1 ' 2 

provides a convenient measure of the polydispersity. 



III. RESULTS 



We now demonstrate the results of the model in two numerical examples, representing two amphiphiles of differing 
hydrophobicity. The parameters of the amphiphiles are given in Table QJ amphiphile B being the more hydrophobic 
of the two. We use here the same two examples whose equilibrium properties have been thoroughly analyzed in Ref. 
[lfj| . Table [J lists for these examples the volume- fraction bounds of the premicellar regime, C2 and C3, along with the 
aggregation numbers at these points, as obtained from the equilibrium theory. 

In Table HT1 we give the values of the free-energy barrier for aggregate dissociation at the lower and upper bounds 
of the premicellar regime, as calculated from Eqs. (HJ) — (I3J) using the parameters of Table [TJ At <f> = Ci the barrier is 
negligible, of order k^T, yet, as <f> increases through the premicellar regime, it becomes much larger than k^T. The 
resulting lifetimes, as calculated from Eq. ([M]) , are given in Table HT1 As an estimate for the molecular time scale we 
have used for both amphiphiles To = 10 ns. (This is the diffusion time of a molecule, having a diffusion coefficient 
of 10 -6 cm 2 /s, along a distance of 1 nm.) Corresponding to the increase in the free-energy barrier, the aggregate 
lifetime increases from milliseconds at the lower bound of the premicellar region to practically indefinite time. As 
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TABLE I: Parameters and equilibrium properties of exemplary amphiphiles. n — number of groups in hydrocarbon tail; wo, 
a, k — parameters of u{m), the free energy of transfer in units of k B T [Eq. (J2J]. £2, C3 — volume-fraction bounds of the 
premicellar regime; m mln (c2), m mln (c3) — aggregation numbers at these boundaries jlOl ] . 



amphiphile n uo 


a K C2 C3 


m min (c 2 ) 


m min (c 3 ) 


A 13 10 
B 20 14 


11 0.08 8.0 x 10~ 4 2.2 x 10 _a 
14 0.05 1.6 x 10~ 5 6.7 x 10~ 5 


53 
118 


60 
128 


TABLE II: Properties of premicellar 
relative width of size distribution; C2, 
used for the molecular time scale. 


aggregates. Fb — free-energy barrier for dissociation; 
C3 — lower and upper bounds of the premicellar regime. 


r m — aggregate 
A value of to — 


lifetime; w — 
10 ns has been 


amphiphile Ft, (02) 
k B T 


-Fb(C3) T m (C2) T m (c3) 

k B T s s 


w{C2) 


Ul(C3) 


A 1.3 
B 0.5 


30.4 2.0 x 10" 4 1.3 x 10 9 
112.5 4.3 x 10" 2 1.8 x 10 43 


0.18 
0.11 


0.15 
0.10 



already noted in Sec. Ill B( our lifetime analysis is strictly valid only for C2 <C <fi <C C3, and, hence, these values should 
be regarded merely as rough estimates. 

The premicellar aggregate lifetime for amphiphile A, scaled by the molecular time To, is depicted as a function 
of amphiphile volume fraction in Fig. [TJ The roughly exponential increase of lifetime with concentration stems 
from the exponential dependence of r m on the barrier height [Eq. (fL4")) ], which is the main source of concentration 
dependence. Two additional contributions to the dependence of r m on </> are included in the prefactor of Eq. (|14p . The 
first, (1 + i 7 "'(m max )/2)~ 1 , comes from the aggregation-number diffusion coefficient, D(m). This factor is practically 
concentration- independent, since in our examples the curvature of the saddle point is small, |F"(m max )| < 0.1, and 
thus D(m max ) ~ (2r ) _1 = const. The second pre-exponential factor in Eq. HI]), \F' '(m min )F"(TO max )| -1 / 2 , depends 
on concentration primarily through \F' '(TO max )|, which is an increasing function of <fi. This factor causes the slightly 
weaker increase of lifetime with <j) at small concentration (Fig. [I]). 




FIG. 1: Lifetime of premicellar aggregates of amphiphile A as a function of amphiphile volume fraction. The lifetime is scaled 
by the molecular time scale to (of typical order of 10 ns), and the volume fraction by C2, the onset of premicellar aggregation. 
The volume fraction corresponding to the cmc (C3) is indicated by an arrow. Parameters of amphiphile A are given in Table U 

In Fig. [1] we see that in the case of amphiphile A, assuming tq ~ 10 ns, the aggregate lifetime reaches the order of 1 
s for <j) ~ 1.4c2, whereas the cmc is at C3 ~ 2.75C2. Thus, the premicellar aggregates remain stable for a macroscopic 
time over a significant part of the premicellar regime. In the case of amphiphile B we find r m ~ 1 s for cf> ~ 1.05C2 while 
C3 ~ 4.2c2, i.e., the aggregates are kinetically stable over a much larger portion (practically all) of the premicellar 
concentration range, as expected for a more hydrophobic amphiphile. 

Figure [5] shows the mean-square fluctuation of the aggregation number for amphiphile A as a function of volume 
fraction. The corresponding relative width of the aggregate size distribution is presented in the inset. The polydisper- 
sity weakly decreases with concentration, i.e., the premicellar aggregates are nearly as monodisperse as the micelles 
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above the cmc. In Table [TT1 we see that the same conclusions hold for amphiphile B. The small polydispersity (around 
10%), as well as the slightly increased value for the less hydrophobic amphiphile (A) , are in agreement with the well 
known trends for spherical micelles above the cmc, as established experimentally |Tlj | and theoretically Q. 
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FIG. 2: Mean-square fluctuation of aggregation number for amphiphile A as a function of amphiphile volume fraction. The 
volume fraction is scaled by C2, the onset of premicellar aggregation. The cmc (03) is indicated by an arrow. The inset shows 
the relative width of the aggregate size distribution, w = (Sm 2 ) 1 '' 2 / (m). Parameters of amphiphile A are given in TablelU 



IV. DISCUSSION 



It has been shown in Sec. IIIII that considerations of aggregate lifetime can reduce the concentration range in 
which premicellar aggregates may be experimentally observable and technologically relevant, compared to the range 
determined from equilibrium considerations alone [lOJ. In other words, the apparent concentration, above which an 
appreciable amount of metastable micelles appears, may be higher than ci- We have demonstrated, nonetheless, 
that kinetic stability (i.e., macroscopic lifetime) still exists in most of the premicellar region. The more hydrophobic 
the surfactant, the wider the range of stability is. These conclusions are in line with results presented in Ref. [2C| . 
Although that study does not deal with premicellar aggregation, it has shown that the dissociation time of micelles 
remains very large even below the cmc. 

Our Kramers-like approach, as already mentioned in Sec. |IIJ relies on two assumptions, which are violated near 
the edges of the premicellar region. The first assumption, of a high free-energy barrier between the metastable and 
stable states, is valid in almost the entire region except close to the lower edge, C2, where the barrier may become of 
order k^T only. (See Table HTl) The resulting short lifetimes, though not accurately accounted for by the theory, are 
of little interest. The second assumption, of a large free-energy difference between the two states, holds in nearly the 
entire range as well, except very close to the upper edge, C3, where, by definition, the frcc-cncrgy difference vanishes. 
The free-energy difference, in units of fepT, becomes large quickly as <fi gets smaller than C3, since the considered 
mesoscopic subsystem of volume V\ contains a large number of molecules (mostly monomers). In addition, correction 
of the theory near C3 by considering a probability backflow from the monomeric to the aggregated state will only 
increase the stability of the latter. Therefore, the deficiencies of the theory at the edges of the premicellar region do 
not affect our main results. 

It should be borne in mind, however, that the stability of premicellar aggregates imply also that high nucleation 
barriers may need to be overcome in order for them to form in the first place. Correspondingly, the more hydrophobic 
the surfactant, the higher these barriers are. (The issue of high nucleation barriers for micelle formation above the 
cmc has been underlined also in Ref. (20|.) Hence, since both the preceding work [l(| and the current one have 
assumed full equilibration, their applicability within reasonable time scales might require in practice either reduction 
of nucleation barriers through heterogeneous nucleation or overcoming them by external means (e.g., agitation or 
sonication) . 

Finally, we have found narrow size distributions of premicellar aggregates, i.e., micelles below the cmc should 
be only slightly more polydisperse than their counterparts above the cmc. (See Fig. [5]) This agrees with the 
monodispersity observed in experiment [6|. Thus, polydispersity does not pose a problem for the applicability of 
premicellar aggregation. 
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